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Abstract 

We study the single transverse-spin asymmetry for the inclusive open-charm production 
in the pp-collision, p^p — > DX, induced by the three-gluon correlation functions in the 
polarized nucleon. We derive the corresponding twist-3 cross section formula in the leading 
order with respect to the QCD coupling constant. As in the case of the semi- inclusive deep 
inelastic scattering, ep^ — > eDX, our result differs from the previous result in the literature. 
We also derive a "master formula" which expresses the twist-3 cross section in terms of the 
gg — > cc hard scattering cross section. We present a model calculation of the asymmetry 
at the RHIC energy, demonstrating the sensitivity of the asymmetry on the form of the 
three-gluon correlation functions. 
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1 Introduction 



Open charm production in inclusive hard processes, such as semi-inclusive deep inelastic 
scattering (SIDIS), ep — > eDX, and the D-meson production in pp collision, pp — > DX, 
is an ideal tool to investigate the gluon distributions in the nucleon. Similarly the single 
spin asymmetry (SSA) in those processes plays a crucial role to reveal the multi-gluon 
correlations in the nucleon [1, 2, 3, 4, 5, 6]. When the transverse momentum of the final 
Z)-meson can be regarded as hard (P T Aq C d), one can analyze the processes in the 
framework of the collinear factorization [7, 8, 9]. In this framework, SSA appears as a twist- 
3 observable and can be represented in terms of the multi-parton correlation functions. 
The purely gluonic correlations responsible for SSAs in the open charm production are 
represented by the "three-gluon" correlation functions. 1 For the quark-gluon correlation 
functions in the nucleon, there have been many studies in the literature in connection with 
SSAs for the light hadron productions [7]- [23], and our understanding on the mechanism of 
SSA has made a great progress. 

In our recent paper [4] we studied the contribution from the three-gluon correlation 
functions to SIDIS, ep^ — > eDX. In that study, we have identified the complete set of the 
three-gluon correlation functions and established the formalism for calculating the twist-3 
single-spin-dependent cross section induced by those functions. The gauge invariance and 
the factorization property of the cross section have been shown explicitly in the leading 
order with respect to the QCD coupling constant. The result of that study differs from 
the previous study in the literature [2] , and we clarified the origin of the discrepancy. In 
our another paper [6], we have developed a novel "master formula" for the three-gluon 
contribution to ep^ — > eDX which expresses the corresponding twist-3 cross section in terms 
of the twist-2 7*<? — > cc scattering cross section, extending the similar formula known for 
the soft-gluon-pole (SGP) contribution to SSA from the quark-gluon correlation functions 
in the nucleon [16, 17]. This formula simplifies the actual calculation of the twist-3 cross 
section and makes its structure transparent, and may be useful to include higher order 
corrections to the asymmetry. 

The purpose of this paper is to extend these studies to the contribution of the three- 
gluon correlation functions to the SSA in the pp collision, 

p^p,S ± )+p(p')^D(P h )+X, (1) 

where S± is the transverse spin vector of the polarized nucleon and P^ is the momentum of 
the D-meson satisfying P% = ml with the D-meson mass m^. The initial nucleons' momenta 
p and p' are in the collinear configuration and can be regarded as massless (p 2 = p' 2 = 0) 
in the twist-3 accuracy. We will derive a corresponding formula for the twist-3 single-spin- 
dependent cross section in the leading order with respect to the QCD coupling constant. 
We will further derive a master formula for (1) which connects the twist-3 cross section to 
the cross section for the gg — > cc scattering. We will also present a model calculation for 

1 ln the framework of the transverse- momentum-dependent factorization, which is useful to describe 
the low- Pt hadron production, the corresponding gluonic effect is represented as the fcj_-dependent gluon 
distribution functions [24, 25]. 
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the asymmetry A N = Aa/a = (cr t — cr^)/(a^ + a^), where represents the cross section 
for (1) with the initial nucleon polarized along S± (—S±), and will obtain a constraint on 
the three-gluon correlation functions, using the recent RHIC data on A^ for the D-meson 
production [26]. 

The remainder of this paper is organized as follows: In section 2, we recall the complete 
set of the three-gluon correlation functions in the transversely polarized nucleon which are 
relevant for our study. In section 3, we derive the twist- 2 unpolarized cross section for the 
process pp — > DX induced by the gluon-density in the nucleon. In section 4, we derive the 
twist-3 single-spin-dependent cross section induced by the three-gluon correlation functions 
for the process (1), applying the formalism in [4]. In section 5, we derive the master formula 
which expresses the twist-3 cross section for (1) induced by the three-gluon correlation 
functions in terms of the cross section for gg — > cc scattering in the twist-2 level. In section 
6, we present a model calculation of An for the D-meson production at the RHIC energy, 
and demonstrate a sensitivity of the asymmetry on the form of the three-gluon correlation 
functions. Section 7 is devoted to a brief summary. In the appendix, we discuss some 
technical aspects in the derivation of the contribution from the initial-state-interaction 
diagrams to the twist-3 cross section. 

2 Three-gluon correlation functions in the transversely 
polarized nucleon 

The twist-3 three-gluon correlation functions in the transversely polarized nucleon were 
first introduced in [1]. Then, as was clarified in [27, 28, 4], there are two independent 
three-gluon correlation functions due to the difference in the contraction of color indices. 
Following the notation in [4], we call those functions 0(xi,X2) and N(xi,x 2 ), which are 
defined from the lightcone correlation functions of the three field-strengths of the gluon in 
the polarized nucleon as 

O a ^( Xl ,x 2 ) = -g(i} 3 J £ J ^e^ 1 e ^ 2 -^)(^|4 C aFr(0)F c ^( / in)Fr(An)b5) 

= 2iM N [0(x 1 ,x 2 )g af> e inimS± + 0(x 2 ,x 2 - x 1 )g fh c? pnS± + 0{x 1 ,x 1 - x 2 )g' ra e PpnS± ] , (2) 
N a ^( Xl ,x 2 ) = -g(i) 3 f ^ J ^ e **i e *(«-* ^(pS\if bca F b ^(0)FJ n (i,n)F: n (\n)\pS) 

'Nixux^g^e^ 1 - - N(x 2 ,x 2 - Xl )g^e apTlSl - - N(x uXl - x 2 )g ia e PpnS A .(3) 



= 2iM 



N 



where = d a A% - A a a + gf abc A%A(! is the gluon's field strength, and we used the 
notation F" n = F^np and e apnS± = e a ^ uX p fl n u S±x with the convention e i23 = 1- d bca and 
jbca are S y m metric and anti-symmetric structure constants of the color SU(3) group, 
and we have suppressed the gauge-link operators which ensure the gauge invariance. p is 
the nucleon momentum, and S± is the transverse spin vector of the nucleon normalized 
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as Sj_ — — 1. In the twist-3 accuracy, p can be regarded as lightlike (p 2 = 0), and n is 
another lightlike vector satisfying p ■ n = 1. To be specific, we set p^ = (p + , 0, 0j_), = 
(0, n~,0j_), and = (0,0, S±). The nucleon mass M N is introduced to define 0(xi,x 2 ) 
and N(xi,X2) dimensionless. The decomposition (2) and (3) takes into account all the 
constraints from hermiticity, invariance under the parity- and time-reversal transformations 
and the permutation symmetry among the participating three gluon-fields. The functions 
0(xi,X2) and N(xi,X2) are real and have the following symmetry properties, 

0(x!,x 2 ) =0(x 2 ,x 1 ), 0(x 1 ,x 2 ) = 0(-X!, -x 2 ), (4) 

N(x 1 ,x 2 ) = N(x 2 ,x 1 ), N(x 1 ,x 2 ) = -N(-x 1 ,-x 2 ). (5) 



3 Twist-2 unpolarized cross section for pp — > DX from 
gluon-fusion 

We first recall the twist-2 unpolarized cross section for the process (1) which is the denom- 
inator of An. It receives main contribution from the cc-creation due to the gluon-fusion 
with the subsequent fragmentation of the c (or c) quark into the D (or D) meson (Fig. 1). 
The corresponding unpolarized cross section can be written as 

h f=c,c 

where S = (p+p') 2 is the center-of-mass energy squared and a s = g 2 / (4ir) is the strong cou- 
pling constant. The c — )■ D (or c — > D) fragmentation function Df(z) and the unpolarized 
gluon distribution in the nucleon G(x) are, respectively, defined as 

£ 7f I ^ e "^( l^(0)|^(^)^>( J D(^)^|^(A^)|0) = (j e + m c ).. + • • • , (7) 

X 

1 1 ^(P\ F r(U)Fr(*n)\p) = -\G{x)sf? + • • • , (8) 

where N = 3 is the number of colors, p c is the momentum of the c (or c) quark fragmenting 
into the D (D)-meson with p 2 = m 2 and w is another lightlike vector of 0(l/p + ) satisfying 
Ph ■ w = 1. Thus p c is related to Ph as = + {m 2 z — m 2 /z)/2w^. is defined 

as g^ = 5^ — p^n" — p"n M . The symbol • • • denotes higher-twist contributions which are 
irrelevant here. The partonic hard cross section ^ g ^ c can be obtained from the 9 diagrams 
shown in Fig. 2 in the leading order (LO) with respect to the QCD coupling constant, and 
is given by 

= in (is " c;p) ( j + * + 4m?s " -it) • (9) 
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Figure 1: Generic diagrams for the twist-2 cross section for pp — > DX induced by the gluon 
densities in the initial nucleons. 




Figure 2: Leading order diagrams for the twist-2 unpolarized hard cross section for <3"^,^ c 
appearing in (6). 

where the invariants for gg — > cc scattering are defined as 

s = (xp + x'p') 2 i = (p c — xp) 2 — m 2 c u = (p c — x'p') 2 — ml, (10) 

and Cf = 

4 Twist-3 cross section for p*p — > DX induced by the 
three-gluon correlation functions 

The twist-3 single-spin-dependent cross section for p^p — > DX induced by the three-gluon 
correlation functions can be obtained by applying the formalism developed for ep^ — >■ 
eDX [4]. The twist-3 cross section occurs from the diagrams of the type shown in Fig. 
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(a) 



(b) 



Figure 3: Generic diagrams which give the twist-3 cross section for p'p — > DX induced 
by the purely gluonic effect in the polarized nucleon (lower blob) convoluted with the 
unpolarized gluon density (upper blob) and the twist-2 fragmentation function for the D- 
meson (middle blob). A pair of circles in each figure represent gg — > cc hard scattering 
amplitudes. Figures (a) and (b) represent, respectively, the contributions from the final- 
state-interaction (FSI) and the initial-state- interaction (ISI). The mirror diagrams also 
contribute. 

3, where the extra coherent gluon is exchanged between the hard scattering part and the 
nucleon matrix element. In Fig. 3, the gluon density G(x') in the unpolarized cross section 
(upper blob) and the fragmentation function Df(z) for the D-meson (middle blob) are al- 
ready factorized. Owing to the symmetry property of the correlation functions (lower blobs 
of Figs. 3(a) and (b)) defined by 

M$?(h, k 2 )=gj J d 4 V e^e^- k ^(pS\A^A x MAm\pS), (11) 

the SSA occurs only from a pole part of an internal propagator in the corresponding hard 
part 

S; b u c x(kuk 2 ,x'p\p c ) J (12) 

where a, b, c are color indices and the momenta k\ and k 2 are assigned as shown in Fig. 
3. Here and below we employ the convention that the QCD coupling constant g associated 
with the attachment of the coherent gluon into the hard part is included in the matrix 
element (11) consistently with the definition of the three-gluon correlation functions in (2) 
and (3). The hard part of the diagrams in Fig. 3 gives rise to the pole contributions at 
X\ = x 2 . (See discussions below.) For those contributions, following the same step as [4] 
in the collinear expansion to S®l c x (ki,k 2 ,x'p',p c ), we eventually end up with the following 
expression for the LO twist-3 cross section induced by the three-gluon correlation function 
(see Appendix): 
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X 



dS%(k u k2,x!lt,PcW 



dk% 



k% — %iP_ 



pole 



^>%M^ c (^i^ 2 ), (13) 



where u^ a = g fJ, a — p^ria, and M^ c (xi,x 2 ) is the lightcone correlation function of the 
field-strengths defined as 



Nd, 



(N 2 



(14) 



with O al3 ' y (xi, x 2 ) and N al3 " t (xi, x 2 ) defined in (2) and (3), respectively. The symbol [• • -] pole 
indicates the pole contribution is to be taken from the hard part. We emphasize that 
even though the analysis of Fig. 3 starts with the gauge-noninvariant correlation function 
(11) and the corresponding hard part S^ c x (ki, k 2 , x'p',p c ), gauge-noninvariant contributions 
appearing in the collinear expansion either vanish or cancel and the total surviving twist-3 
contribution to the single-spin-dependent cross section can be expressed as in (13), using 
the gauge-invariant correlation functions (2) and (3). 

- pole 



The pole contribution to the hard part 



dsp x (k u k 2 ,W,Pc)p/dkz 



=x,p 



occurs 



from two types of diagrams shown in Figs. 3 (a) and 3(b), which are refered to as the 
final state interaction (FSI) and the initial state interaction (ISI), respectively, by the 
parton lines to which the coherent gluon is attached. Figs. 4 (a) and 4(b) show the LO 
diagrams contributing to the hard part S°^ x (ki, k 2 , x'p',p c ) in Figs. 3(a) (FSI) and 3(b) 
(ISI) , respectively. The mirror diagrams of Fig. 4 also contribute. The poles are produced 
from the bared propagator, and gives rise to the (^-function at x\ = x 2 in the collinear 
limit (ki — > Xip), hence the poles are refered to as the soft-gluon-pole (SGP). Other pole 
contributions cancel among each other by taking the sum of the whole diagrams. 



By calculating dS^ c x (ki, k 2 , x'p' ,p c )p x / dk^ 

coefficient tensors in the decomposition of (2) and (3~ 
dependent cross section as [5] 



pole 

from Fig. 4 contracted with the 
, one obtains the twist-3 single-spin- 



dAa_ 

h cPP h 



a s M N ir PhpnSj 
S 



E 

f=cc 



dr' 

d -^G(x') 



dz „ , . 



dx . ~ . 1 
— 5(s + t + u) — 

X zu 



2Q(x,x) 



x 



dx x x x 



+ 



d_ 

dx 



N(x,x) - 



2N(x,x) \ &m + (d_ N , Xj Q) _ 2iV(x,0) \ a N2 + N(x,x) „ m + N(x,0) ^ N4 
x J \dx x J x x 



(15) 
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(a) 




Figure 4: The LO diagrams for the partonic hard part S^ c x (ki, k 2 , x'p' , p c ) for the twist-3 
cross section. Diagrams in (a) represent the FSI contribution and those in (b) represent 
the ISI contribution. 



where S- — 1 and 5 P 
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■1. The partonic hard cross sections are given by 
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(16) 
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and 

\ m ( 1 t 2 + u 2 1 u 1 s \(r> „ 2 A 2 _ 4m 4 J 2 \ 

s« = _ (^LtS + -LA _ _L_ ') (p + ^ + Smls _ , 

\C F s 2 tu C F st 2 N 2 C F t 2 u \ c tu 

< (17) 

/ 1 t 2 + U 2 11 1 5 2 \ / 4 _ „ 2 __x 

* m =~ few + C^P - N^M ^ - ^ • 

The hard cross sections associated with the first term in the first parentheses in (16) and 
(17) come from ISI (Fig. 4(b)), and those associated with the second and the third terms 
in the same parentheses come from FSI (Fig. 4(a)). As in the case of ep^ — > eDX, the cross 
section in (15) receives the contribution from the four functions 0(x,x), O(x,0), N(x,x) 
and N(x, 0). Unlike the case of SIDIS, presence of ISI gives rise to the different hard cross 
sections for O and N functions. From (15), it is clear that the process p^p — >■ DX itself is not 
sufficient for the complete separation of the four functions. For the separation, the process 
ep^ — > eDX serves greatly, since it has five structure functions with different dependences 
on the azimuthal angles to which the four functions contribute differently [4] . For the 
massless quark fragmenting into a light hadron (i.e. m c — > ), one has o- 03 . 04 ^ 3 ^ 4 _ > o, 
a° l = o° 2 and a N1 = —a m . Therefore the three-gluon correlation functions appear in the 
combination of x(d/dx)(0(x,x) + O(x,0)) — 2(0(x,x) + 0(x, 0)) and x(d/dx)(N(x, x) — 
N(x, 0)) - 2(N(x, x) - N(x, 0)) at m c = 0. 

Our result in (15) differs from a previous work [3]: The result in [3] is obtained from 
(15) by omitting the terms with a 02 ' 04 and cr N2 > N4 and by the replacement 0(x,x) — > 
0(x, x) + 0(x, 0) and N(x, x) — > N(x, x) — N(x, 0). This difference originates from an ad- 
hoc assumption in the factorization formula in [3, 2]. We emphasize the appearance of the 
four different contributions with {0(x, x), 0(x, 0), N(x, x), N(x, 0)} is a consequence of the 
symmetry property implied in the decomposition (2) and (3), in particular, the different 
coefficient tensors in front of 0(x, x) and 0(x, 0) (likewise for N(x,x) and N(x,0)) at 
x 1 = x 2 = x lead to different hard cross sections for the above four functions. See [4] for 
more details. 



5 Master formula for the three-gluon contribution to 

p^p -> DX 

5.1 Connection between the twist-3 cross section and the gg — > cc 
scattering 

To obtain the twist-3 cross section based on (13), one has to calculate the derivative of the 
hard part [dS^. c x (ki, k 2: x'p' ,p c )p x ldk 2 ] pole from Fig. 4 contracted with the coefficient 

^ ki=Xip 
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tensors in the decomposition of (2) and (3). This calculation produces lots of terms at the 
intermediate step and is extremely complicated. Alternatively, application of the "master 
formula" developed for the contribution of the quark-gluon correlation functions [16, 17] and 
also for the three-gluon correlations for ep^ — > eDX [6] provides us with a more transparent 
and simpler method to calculate the cross section. To extend the method to the contribution 
of the three-gluon correlation functions for p^p — > DX, we first note that the diagrams in 
Fig. 4 are obtained by attaching the extra gluon-line to the external-lines of the twist-2 
hard part in Fig. 2, and the pole contribution is given by the propagator next to the 
vertex to which this extra-gluon line attaches. Because of this structure, the derivative 
can be performed by keeping the structure of the hard part corresponding to those in Fig. 
2 almost intact. Based on this observation, one can obtain the master formula also for 
the three-gluon contribution. To be specific, we consider the case in which the c-quark 
fragments into the .D-meson below. 

To present the result, we first define the hard part for the unpolarized cross section 
shown in Fig. 2 as 7i^ b (xp,x'p',p c ), where \iv and ab are, respectively, the Lorentz and 
the color indices for the gluon line with the momentum xp. Those indices for the gluon 
line with the momentum x'p' are already contracted to factorize G(x') in (6). With this 
convention the partonic hard cross section (J gg ^ c in (6) is related to "H^ ab as 



[s, t, u, m 2 c )5 (s + t + uj = 



(N 2 - 1) 



■^f)^ ab (^y,p c ). (is) 



As is shown below the hard part for the twist-3 cross section has a simple relation with this 

In [6], we have shown that the twist-3 hard cross section for ep^ — > eDX induced 
by the three-gluon correlation functions can be obtained from the Born cross section for 
the 7*g — > cc scattering. There the SGP contribution occurs from the FSI. Accordingly, 
the FSI contribution for p^p — > DX shown in Fig. 4(a) can also be expressed in terms 
of the Born cross section for the gg — > cc scattering. We write the FSI contribution to 
S^l° x (ki, k 2 , x'p',p c )p x in (13) as S^ c (ki, k 2 , x'p',p c )p x . Then, following the same procedure 
as [6], one can show that the hard part for the FSI is given by 



dS^ a x C (k 1} k 2 ,x'p' \p c )p } 



dk% 



pole 



1 



-i pole 



x i — x 2 + ie 
1 



xi — x 2 + ie 



d p ca p x d 



\ dp a c p ■ p c dp 

pole d 



Kf c {x lP ,x'p',p c ) 



^n^ix^x'p'^x 



(19) 



where 'H^ bc (xip,x'p',p c ) is obtained from H^^^Xip, x'p',p c ) simply by adding the extra 
color matrix t c in the same place where the coherent gluon line is attached in Fig. 4(a). 



10 



Here one needs to be cautious in taking derivative with respect to p°: In the expression 
after the first equality of (19), the on-shell limit p 2 = m 2 should be taken after performing 
the derivative with respect to p a c . For the derivative in the expression after the second 

equality of (19), the form p^ = (p+ = mc ^- c± , p~ , p C A^j should be used for p c , i.e., on-shell 

condition for p c should be used by regarding p+ as a dependent variable of p~ and p c ±. 

One can also derive the similar relation for the ISI contribution. We write the ISI 
contribution to S"^ c x (ki, k 2 , x'p', p c )p x in (13) as S I '^ c (k\, k 2 , x'p', p c )p x - For the ISI diagrams 
in Fig. 4(b), the coherent gluon couples to the initial gluon-line of the diagrams in Fig. 2 
through the three-gluon coupling. One can still apply the same method as [6] , and obtains 
for the ISI contribution as 



dS I ^ c (k u k 2 ,x'p',p c )p* 



dk% 



hi — %iP 



pole 



x 2 — x i + 



poic ( Q p' p x d \ 

\d(xW) ~~ p a -p'd(x'p^) n ^ bC{XlP,X ' P '' Pc 



-1 



x 2 — x\ + ie 



pole 



U^ix^x'p'^), 



d(x'p' ( 



(20) 



where 'H I ^ bc (x 1 p,x'p' ',p c ) differs from H^ b (xip,x'p',p c ) only with its extra color index c 
associated with the attachment of the coherent gluon line in Fig. 4(b). As in (19) the 
on-shell limit p' 2 — > should be taken after carrying out the derivative in the expression 

after the first equality in (20), and the on-shell form p /fl = ^p' + = ^h,p'~ ,p'±^j should be 

used in the expression after the second equality of (20). In (20), we first make p"]_ ^ 
in taking the derivative and then take the p"]_ — > limit to consider the cross section in 
the frame where p and p' are collinear. Inserting (19) and (20) into (13), one obtains the 
single-spin-dependent cross section as 



x 



d 



d 



n^ bc (xp,x'p', Pc ) 



(21) 



The hard part "H^ a6c and "H^* c contain the factor 5 ((xp + x'p' — p c ) 2 — m 2 ) — 5 (s + t + u 
as an on-shell condition for the final unobserved c-quark. For convenience we separate this 
5-function and introduce the following functions by taking the color contraction: 

(AT2 _^(AT2 _ ^ n a'p bC (xP, X'p', Pc) = (X P , x'p', Pc )5 (§ + t + u) , 

AT(^-l) ^ C(m X ' V ' ' Pc) = H ™ {XP > X ' V ' ' Pc)5 ( S + * + fi ) ' 
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{N 2 - N 1){N 2 - x V,Pc) = H { a y\xp, x'p',p c )5 (s + t + u), 



Using these forms in (21), one can write the cross section as 
dA(x [ dx', f dz f dx . , 



(22) 



X 



0? 7 (x, x) I -^H { a F /\xp, x'p',p c ) - j H^ d \xp, x'p',p c ) \S (s + t + uj 



d 



dp, 

d_ 

dp, 



d(x'p^) 



+Nf (x, x) \ -^H%^ (xp, x'p',p c ) - ^A-^/) (xp, x'p', p c ) | 5 (§ + t + uj 



,(23) 



where O^^x, x) and N'^ 1 (x,x) are the functions obtained by setting X\ = x 2 = x in (2) 
and (3): 



O a /\x, x) = 2iM N 0{x, x)gl p e^ n ^ + 0(x, 0)(^ 7 e Qpri ^ + g^e 



at/3 jypnS± 



^l £ apnS ± , 7" c /3pri5 ± ^ 



Nf 1 (x,x) = 2iM N \N(x, x) gfe lpnSl - - N(x, 0) (g^e apnSl - + gl a e^ nS± )} . (24) 



We remind the derivatives d/dpj and d/d(x'p' 1 ) in (23) also hit S [s + t + uj . The relation 
(23) shows that the partonic hard cross section for 0(x, x) and 0(x, 0), in general, differ 
from each other, and likewise for N(x,x) and jV(x, 0). 



5.2 Contribution from 0(x, x) and N(x,x) 

We first consider the contribution in (23) occuring from 0(x,x) and N(x,x) in (24). Using 
the functions (22), we write the corresponding hard part as 

H^ j \xp,x'p',p c )gfe^ = K^\s,t,u,m 2 c )e^ s \ 

H { Jf (xp, x'p', p c ) gf e ^ nS ^ = K {I ' j) (s, t, u, m 2 c )e lpnS ^ , (25) 

for j = d, f, where we have used the fact that the scalar functions K^ F ^ and K^ 1 ^ (j = 
d, f) become the functions of s, t, u and m 2 c . For the scalar functions K^'fi and 

5 (s + 1 + uj , one can perform the derivative with respect to pj and x'p' 1 in (23) through 
that with respect to u as 

^K^\s, t, u, m 2 c )5 (§ + t + u) = -2p n (0 ^-K^\s, t, u, m 2 c )5 (s + t + u) ,(26) 
K^\s, t, u, m 2 c )5 (s + t + u) = -2p n -^K^\~s, t, u, m 2 c )5 (s + t + u), (27) 
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for j = d, f, where we have used the fact that p+ and p' + are the dependent variables (as 
noted after (19) and (20)), and have set p'\ —¥ after taking the derivative. 2 Using (26), 
one obtains for the FSI contribution with 0(x,x) in (23) as 



d 



dp, 



0^ypnS± 



= _ 2e p cP nS ± J ^0(x,x) (Vj ^K^ d \S,t,u,m 2 c )5(s + t + u) 



-2e 



p c pnS± 



dx /S\ fdK™ xdK (F ^ 2K^\ . K (F ^ dO(x,x) 

+ -— ^ = — }0(x,x) + ^ Z7 —x- ' 



x 2 \tJ [\ du u dx 
x5 (s + t + u^j 



u 



u 



dx 



dx f s 
'x 2 If 



fdK™ tdK™ sdK™ 
\ du 



+ 



+ 



_2 e PcP nS ± 



dx f s 
x 2 \tu 



—m„ 



u 
dm 2 . 



u dt u ds 

5(§ + t + u) 



0(x, x) 



dO(x,x 

H z — x 



dx 
0(x, x) + 



R{m L d °M -2Q(x,x) 



x 



5 (s + i 



u 



(28) 



In the first equality of (28), we transformed the derivative hitting 5 (§ + 1 + uj into the 
derivative with respect to x and performed the partial integration. In the last equality we 
have used the relation 



'_9 ~d „d 
s— + t— - + u-— + m 



d 



ds dt 



du 



'dm 2 



K^\s,t,u,m 2 c ) = 0, 



(29) 



for B = F,I and j = d, f, resulting from the scale-invariance property for the dimensionless 
function K^ B ' j \s,t,u,m 2 c ) = K^ B ' j \\s, \t, \u, \m 2 c ). Similarly to (28), by using (27), one 
obtains for the ISI contribution with 0(x, x) in (23) as 



= 2e p cP ns ± f <^ x )d K {Wf S ,t,u,m 2 )5(s + t + u) 
J x 2 du v ' 



2 For the FSI, one can set p'\ = from the beginning. 
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= 2e PcPnS ± f 

J X 2 U 



-m c 



2 dK^ 
c dm 2 



0(x,x) + K( I 'V 



x 



dO(x, x) 



dx 



20(x,x)\ S(§ + t + uj. 



(30) 



By the replacement K^ B ' d \s,t,u,ml) — > ,t,u,m 2 ) (B = F, I) and 0(x,x) — > 

N(x,x) in (28) and (30), one obtains the formula for the N(x,x) contributions. From (28) 
and (30), one can make the following important observations for the 0(x,x) and N(x,x) 
contributions: 

(1) The partonic hard cross sections for x-^0(x ) x) and 0(x,x) are connected by the 
simple relation. In particular, in the m c — > limit, they contribute in the form 
x-^0(x,x) — 20(x,x). The same relation holds for x-^N(x,x) and N(x,x). 

(2) The partonic hard cross sections K^ B ^ (B = F, J, j = d, f) defined in (25) are 
obtained from U F f c and U^ bc by the same Lorentz contraction as the unpolarized 
cross section in (18). Therefore the contribution to K^ B ^ from each diagram differs 
from those for &gg^ c only in the color factors. 

These features are the extension of those obtained in [17] for the SGP contribution of 
the quark-gluon correlation function with massless partons to the case of the three-gluon 
correlation functions with massive partons in the final state. 

5.3 Contribution from O(x,0) and N(x,0) 

Next we consider the contribution from 0(x, 0) and N(x, 0) in (23), which arise from the 
second terms in (24). For this purpose, we introduce the two fixed vectors X^ = (0, 1, 0, 0) 
and = (0, 0, 1, 0), and write 



Then the derivative d/dpj hitting the FSI hard part in (23) for 0(x, 0) and N(x, 0) can be 
written as 



Pi = —X^X 1 — Y^Y 1 



(31) 




= -X^H { a F f\xp, x'p',p c )5 (s + t + u) (Ve apri5 ^ + X a e PpnS ±) 

-Y^-^H { a F /\xp, x'p',p c )S (s + t + u) (Y^e apnSl - + r a e* nSl ) . (32) 



To perform the derivatives in this equation, we introduce the scalar functions J 1)2 J (s, t, it) 
(j = d, f) by the decomposition: 



H { a F ' j \xp,x'p',p c ) (xfe apnS± +xV pn ^) 

= 4 F ' j) (s, t, u, m 2 c ) (p c ■ X) t PcPnS ^ + 4 FJ) (s, t, u, m 2 c )e XpnS± , (33) 
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where we used the kinematic condition p]_ = p'\ = 0. One also obtains the similar decom- 
position for H { a F ' j \xp,x'p',p c ) (yPe apnS± + Y a e^ pnS± ^j by the replacement X ->■ Y in (33) 

with the same functions ■ Then the derivative in (32) can be performed with the help 
of (26) as 



d 



j[ F ' l] {§, t, u, m 2 c ){ Pc ■ X)t PcPnS ^5 (s + t + «)] + (X ->• Y) 
= 3e PcPnS ^J (F ' J \s, t, u, m 2 c )5 (s + t + u) 



+2p c \^ pnS - (!) ^ (Jp\~s, t, u, ml)5 (5 + * + «)), (34) 



dp, 



4 Fd \§, t, u, m 2 c )e XpnS ^5 (s + t + u)] + (X -> F) 



s 9 



(35) 



Using these results and the relation p 2 ± = y — m 2 in (32), and following the same procedure 
leading to (28), one obtains the FSI contribution with 0{x, 0) and N(x, 0) in (23) in terms 
of j[ F i j) in (33) as 

dx d 



+ (Hl F > d) ^ H {FJ \ O(x,0) ^ -N(x,0) 



c p c pnS_ 



J 



dx 



x. 



34 F > d) O(x,0) + 2 2 r-^-ml 
tu V s 



{ \ -4 F4) - ™f4^\ 0(x, 0) + J™ L d0 ^) _ 20(x, 0) 



2s 
tu 



dA F ' d) 



-mi^rO(x,0) + J?»\x 



{F4) ( dO(x,0) 



dx 



dx 



-2O(x,0) 



5(s + t + u) 



+ (j$ A ->J$ f) , O(x,0)^-N(x,0 



(36) 



The ISI contribution in (23) with 0(x, 0) and N(x,0) can also be obtained following 
the same procedure as above. Similarly to (33), one can decompose the ISI hard part as 



H^ j \xp,x'p',p c ) (xPe apnS ± +xV pn ^) 
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= j[ IJ \s, t, u, m 2 c ) (p c ■ X) e PcPri5 ^ + 4 IJ \s, t, u, m 2 c )t XpnS ^ 

+4 I ' 3) {S, t, u, m\) (x'p' ■ X) t MmSl - + Jl /J) (§, t, u, ml) (p c ■ X) x'e p ' pnSl - (37) 

for j — d, f, where we ignored the terms which vanish in the limit p n ± — > after taking the 
derivative dj 'dix'p' 1 ). Similar relation can be written down for 

H^ j) (xp,x'p',p c ) (Y l3 e apnS± + Y a e (3pnS± y using the same functions Compared with 

(33), note the existence of the J^ 1 ^ terms in (37), since one has to take the derivative with 

respect to p'\ before taking the p'\ — > limit. With these J12J4, one eventually obtains 
the ISI contribution in (23) with O(x,0) and N(x, 0) as 

dx d 

~~2 



X 



-^^H^ d \xp,x'p\p c )5 (s + t + u) ^ e ^s ±+gTe , pnS ^ 0{x Q) 
+ {H^^Hiy\ O(x,0) ^ -N(x,0) 



p c pnS 1 



dx 



2 tu 2 
-\ — -m c 
u V s 



-J^(x d °[ X ; 0) -2O(x,0) 



dx 



2 

+- 

u 



-m 



2 dJ 2 



(i,d) 



drni 



-O(x,0) + [x 



dO(x,0) 
dx 



2O(x,0) 



-(2jf ' d) + Ji 14) )0{x, 0)] 5 (s + t + u) 
+ (^Sl J$$a, 0(x, 0) -N(x, 0)) . 
5.4 Twist-3 cross section from the gg — > cc scattering 



(38) 



Using the results (28), (30), (36) and (38) in (23), one obtains the final result for the 
total twist-3 single-spin-dependent cross section for p^p — > DX induced by the three-gluon 
correlation functions as 



dx 
1 



d 6 Ph S J x J z 6 J 

'-0(x, x) - (x d °^ X) - 20(x, x)^j j 

-0(s, x) + {x^^- - 20(x, *)) j 



2s 
tu 



dml 



2 f 2 dK^ 
u 



' dm 2 c 
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+ (0(x, x) ->• N(x, x), K^ F,f \ ->■ K^' f ^ 



+34 F ' d) O(x,0) 




x 



j(F,d) 



(F,d) ' 



<9m 2 



(F,d) / dO(x,0) 



eke 



-2O(x,0) 



W«0 („ d ( x i Q ) 
dx 



2O(x,0) 



+ (O(x,0)^-N(x,0), J[f ] j[ F 2 f) ) 




j[ I4) +mf-§^\0{x^)-J[ 



™ (x d -^^-2O(x,0) 



+- <( -m; 



dml 



O(x,0) + Jt > [x 



(i,d) ( dO(x,0) 



dx 



-2O(x,0)\ -(24 



(j,d) 



+ 4 I ' d) )O(x,0) 



+ (O(x,0) ->• -iV(x,0), 



j{I,d) 
J l,2,3,4 



v. 7(/J) 

^ J l, 2,3,4 



S + t + U 



(39) 



where K^\s, t, u, m c 2 ) (5 = F, I, 3 =d,f), J§ j) (5, t, u, m 2 c ) (j = d, f) and j[] 2 %{s, t, u, m c 2 ) 
(j = d, f) are the functions defined, respectively, in (25), (33) and (37) and can be calcu- 
lated from the twist-2 diagrams in Fig. 2. By the direct calculation of these functions, we 
found that j[ F '^ and (j = d, /) are of 0(m 2 c ) and thus vanish in the m c — > limit, and 

that K^ B '^ = J^'^ (B = F, I, j = d, f) at m c = 0, which is consistent with the result in 
(15). For the twist-3 cross section for the D-meson production, the sign of the contribution 
from 0(x ) x) and O(x,0) should be reversed in (39). The result calculated from (39), of 
course, agrees with (15) which were obtained by the direct calculation of Fig. 4. 

The origin of the above master formula (39) is the relations in (19) and (20). Although 
they were derived in the LO QCD, the derivation was based on the quite general structure 
of the diagrams for the SGP contribution at x\ = x 2 shown in Fig. 3, i.e., the coherent- 
gluon line is attached to the external parton line of the twist-2 diagrams in Fig. 1. As long 
as this structure is kept after including higher-order corrections, the relations (19) and (20) 
hold. Therefore we expect that the formula (39) will become a powerful tool to include 
high-order corrections to the twist-3 SSA [6] . 
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6 Numerical calculation of the asymmetry 



As is shown in (15), four nonperturbative functions 0(x,x), N(x,x), O(x,0) and N(x, 0) 
participate in the twist-3 cross section for A®- Unlike twist-2 parton distributions, twist-3 
multiparton correlation functions do not have probability interpretation and thus cannot be 
constrained by a certain positivity bound. 3 They have to be determined by comparing the 
calculated SSAs with experimental data, or by some nonperturbative techniques in QCD. At 
present there is no information on these functions. Preliminary data on A® by the PHENIX 
collaboration [26] suggests |A^| < 5 % in the region \xf\ < 0.1 at y/S = 200 GeV. Since 
the unpolarized cross section for the pion production at RHIC has been well described 
by the next-to-leading order (NLO) calculation in the collinear factorization [30, 31, 32], 
comparison of the asymmetry calculated by our twist-3 cross section with the RHIC data 
will be the first step for determining the magnitude of the three-gluon correlation functions 
in the nucleon. 4 Here we present a simple model calculation of the asymmetry at the 
RHIC energy taking into account of the preliminary RHIC data. 



100 




Figure 5: Unpolarized cross section for pp — > DX by the gluon fusion process at the RHIC 
energy = 200 GeV and P T = 2 GeV. 



To see the relative importance of each term appearing in (15), we assume the same form 
for the four nonperturbative functions as 5 

0(x, x) = 0(x, 0) = N{x, x) = -N(x, 0). (40) 

As a functional form of these functions, we employ the following ansatz: 

Model 1 : 0(x,x) = K G xG(x) } (41) 



3 For a comprehensive review on the positivity bounds of parton distributions, see [29]. 

4 The size of the NLO correction to the twist-3 cross section in (15) could be different from that for the 
twist-2 unpolarized cross section, which may lead to significantly different A®. One thus should take the 
present calculation as only an estimate of the order-of-magnitude. 

5 Minus sign is introduced for N(x,0), since a N2 has an opposite sign compared with jj 01 ' 02 '^ 1 . 
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Model 2 



0(x,x) = K' G y/x~G(x) 



(42) 



where G(x) is the twist-2 unpolarized gluon density, and Kq and K' G are the constants which 
we determine so that the calculated asymmetry is consistent with the RHIC data. Since the 
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Figure 6: (a) Contribution to A® from the 8 components proportional to o 0l,02 > N1 ' N2 (left) 
and 0-03,04,^3,^4 ( r j g ht) in (15) obtained by using the model 1 in (41) with K G = 0.002. In 
the left figure, the non-derivative terms contributing with also plotted by 

thin lines labeled by 0(x, x), 0(x, 0), N(x, x) and N(x, 0). (b) The same as (a) but for the 
model 2 in (42) with K' G = 0.0005. 



three-gluon correlation functions are completely independent from the gluon density, the 
above parametrization is a very crude approximation and the result below should be taken 
only as an estimate of the order-of-magnitude for the three-gluon correlation functions. 
But they are useful to get the shape and magnitude of the three-gluon correlation functions 
relative to the gluon density. Note that the above two models monitor the sensitivity of A 1 ^ 
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to the small-x behavior of the three-gluon correlation functions. For the numerical calcula- 
tion, we use GJR08 distribution [33] for G(x) and KKKS08 fragmentation function [34] for 
Df(z). We also assume the same scale dependence for 0(x,x) etc as G(x) for simplicity. 
We calculate A^ for the D and D mesons at the RHIC energy of VS = 200 GeV and the 
transverse momentum of the D-meson Pt = 2 GeV with the parameter m c = 1.3 GeV by 
setting the scale of all the distribution and fragmentation functions at \i = J P? + m% 

For completeness, we first show in Fig. 5 the unpolarized cross section for pp — > DX 
based on the gluon fusion process (6) at = 200 GeV and P T = 2 GeV. This will 
constitute the denominator of A^, in our calculation below. 




Xp Xp 

(b) 

Figure 7: (a) A 1 ^ for the D° (left) and D° (right) mesons for Model 1 in (41) with Kq = 
0.002. (b) A® for the D° (left) and D° (right) mesons for Model 2 in (42) with K' G = 0.0005. 
Bars denote the RHIC preliminary data taken from [26]. 

Fig. 6 shows the results for the contribution to A® from each term in the twist-3 cross 
section shown in (15) with Kq = 0.002 for the model 1 (Fig. 6(a)) and K' G = 0.0005 for 
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the model 2 (Fig. 6(b)). For both models, one sees that the nonderivative terms accompa- 
nying the hard cross sections o- 03 - 04 ^ 3 ^ 4 [ s negligible compared to the contributions from 
£.oi,02,ni,N2 j n f|g ures f Figs. 6 (a) and (b), we have also plotted the contribution 

from the nonderivative terms proportional to (j 01 >° 2 , N1 , N2 ? which shows that the derivative 
contributions dominate in these terms at large xp (> 0), while the effect of the nonderiva- 
tive term becomes important at xp < and even dominates the asymmetry for the model 
1. One also sees that a 01 and a° 2 give rise to numerically very close asymmetries at the 
RHIC energy for the two models, and likewise for a N1 and a N2 . The asymmetries caused 
by a 01 ' 02 and cr N1,N2 are also similar. 

Fig. 7 shows the result for A® for the D and D mesons including all the contributions 
in (15) together with the preliminary data by the PHENIX collaboration [26]. Because the 
sign of the contribution from {0(x, x), 0(x, 0)} changes between D and D as shown in 
(15), {0(x, x), 0(x, 0)} and {N(x, x), N(x, 0)} contribute to the asymmetry constructively 
(destructively) for the D (D) meson, leading to a large (small) A® for D (D). If one reverses 
the relative sign between O and N from (40), the result for the D and D mesons will be 
interchanged. The values K G = 0.002 and K' G = 0.0005 have been chosen such that the 
calculated asymmetries does not overshoot the data for A N for the D-meson (left figures in 
Figs. 7(a) and (b)) under the assumption (40). By comparing the results for the models 1 
and 2 in Fig. 7, one sees that the behavior of the asymmetry at xp < depends strongly 
on the small- x behavior of the three gluon correlation functions. Therefore A^ at xp < 
is useful to get constraint on the small- x behavior of the three-gluon correlation functions. 

As we saw in the left figures of Fig. 6, we found the relations a 01 ~ a 02 and a N1 ~ a N2 . 
This means that the combinations 0(x, x) + 0(x, 0) and N(x, x) — N(x, 0) can be taken as 
good effective three-gluon correlation functions determining A^'s at RHIC energies. From 
the left figures in Figs. 7 (a) and (b), if the ansatz (41) or (42) is a reasonable assumption 
for the x-dependence of the three-gluon correlation functions, K G = 0.002 or K' G = 0.0005 
can be taken as a modest upper bound corresponding to the assumptions (41) and (42). 
Therefore we may set the upper bound for the combination as 

\0(x,x) + 0(x, 0)| < 0.004 x G(x), \N(x, x) - N(x, 0)| < 0.004a; G(x), (43) 

and 

\0(x,x) + 0(x, 0)| < 0.001 y/xG(x), \N(x,x) - N(x, 0)| < 0.001 y/xG(x), (44) 

although extraction of the separate constraint on the four functions is not possible. 6 We 
remark that even though the RHIC data suggests small A® at \xp\ < 0.1, it can be much 
larger at \xp\ > 0.1 depending on the behavior of the three-gluon correlation functions in 
the large and small x regions as shown in Fig. 7. 

6 In [5], we extracted a stronger constraint for the upper bound of \0(x,x) + O(a;,0)| and \N(x,x) — 
iV(;r,0)|, since we assumed \A®\ < 5 % in the wider region of xp, while the data showing \A®\ < 5 % is 
only in the region \xp\ < 0.1 [26]. Therefore, at present, a weaker upper bound in (43) and (44) is more 
appropriate. 



21 



7 Summary 

In this paper we have studied the SSA for the open-charm production in the pp collision, 
p^p — > DX, based on the twist-3 mechanism in the collinear factorization. Since the three- 
gluon correlation functions in the transversely polarized nucleon play a dominant role in 
giving rise to SSA for this process, we have derived the corresponding twist-3 single-spin- 
dependent cross section in the leading order QCD. As in the case of our previous study 
on ep^ — > eDX, our result differs from the existing result in the literature. We have also 
derived the master formula which shows that the corresponding twist-3 cross section can be 
obtained from the hard part for the gg — > cc scattering in the twist-2 level. The use of this 
formula simplifies the actual calculation and is useful to make the structure of the twist-3 
cross section transparent. We expect that this master formula is useful for the inclusion 
of the higher-order corrections to the cross section. We also presented a model calculation 
of the asymmetry A® in comparison to the preliminary data obtained at RHIC. We have 
shown that at xp < is sensitive to the small-x behavior of the three-gluon correlation 
functions, and have given a modest upper limit on those functions. 
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A Ward identity for the initial state interaction 

To derive the cross section for p^p — y DX, one has to analyze 



where M^ c (k\, k 2 ) is defined in (11) and S°^ c x (ki, k 2 , x'p',p c ) is the corresponding hard part 
as shown in Fig. 3. As was shown in [4], in order to be able to obtain the cross section 
from (13), it is essential that S"^ c x (ki, k 2 , x'p',p c ) satisfies the Ward identities 



kis;l c x (k u k 2 ) = 0, kzs&fa, h) = 0, (k 2 - h^szlKh, k 2 ) = 0. (46) 



It is easy to see that the hard part for the sum of the FSI diagrams shown in Fig. 4(a) 
satisfies (46) due to the on-shell condition for the bared quark-lines. However, the hard 
part for the ISI diagrams in Fig. 4(b) does not satisfy (46). This is because the polarization 
tensor for the gluon line producing the SGP (bared gluon line in Fig. 4(b)) is taken to be 
—g aT in the Feynman gauge, which contains the contribution from unphysical polarizations. 
Nevertheless, one can calculate the hard cross section for the ISI from (13). To show this, 




(45) 
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we define the new hard part for the ISI, S I ^^ c (ki,k 2 ,x'p',p c ), which is obtained from the 

original hard part for ISI, S I ^ x bc (ki, k 2 ,x'p' ,p c ), by replacing the polarization tensor for the 
bared gluon-propagator in Fig. 4(b) as 

~9ar -> -gar H , (47) 

q-p 



where q is the momentum carried by the bared gluon-line. One can show that in the twist-3 

" <I,abc 

d 4 ki f d A k 2 



accuracy this S I '^ c (ki,k2,x'p',p c ) satisfies the relation 



(27r)J j ^yS^ik.M.x'p^p^M^M) 



I 



<l h J^SffikuhMpJM&^kt), (48) 



(27r) 4 7 (2tt) 



where M^(k±, k 2 ) is defined in (11). The extra terms in S^^°(ki, k 2 ,x'p',p c ) compared to 



Suv\ C (kii k 2 ,x'p',p c ) which occur due to the replacement (47) can be shown to vanish with 
the help of the Feynman gauge condition for the matrix element: 

k^M^ih, k 2 ) = 0, k 2u MZ\h, k 2 ) = 0, (k 2 - k^M^ih, k 2 ) = 0. (49) 

From the above relation, one can use S^ c (ki, k 2 ,x'p',p c ) for the calculation of the twist-3 
cross section. Since S I $ c (ki,k 2 ) satisfies 

WS^fih, k 2 ) = 0, k'S^fih, k 2 ) = 0, (k 2 - h^fih, k 2 ) = 0, (50) 

the resulting twist-3 cross section from the ISI takes the form of (13) with S^ x (k\, k 2) x'p 1 ', p c ) 
replaced by Sfflkuh 

,x'p',p c ). Finally, one can show the relation 



9S I ^(k 1 ,k 2 ,x , p / ,p c )p x 



^ a u^u%M^ bc { Xl ,x 2 ) 

ki=Xip 



dk c 2 

_ dS I l £^{k 1 ,h 2 ,x'p',p < )p x 



dk% 



W^&i,^), (51) 



7 

hi = x ip 



where the right-hand-side is calculated with the original hard part for the ISI. This way 
one can obtain the contribution of the three-gluon correlation functions to the twist-3 cross 
section from (13). 
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